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o\ ; 

^ We consider an exact expression for the 6j-symbol for the isosceles tetrahedron, involving SU(2) 

, group integrals, and use it to write the two-point function of 3d gravity on a single tetrahedron as a 

' group integral. The perturbative expansion of this expression can then be performed with respect 

to the geometry of the boundary using a simple saddle-point analysis. We derive the complete 
expansion in inverse powers of the length scale and evaluate explicitly the quantum corrections 
up to second order. Finally, we use the same method to provide the complete expansion of the 
isosceles 6j-symbol with the explicit phases at all orders and the next-to-leading correction to the 
Ponzano-Regge asymptotics. 



(N 



O 



X 



Introduction 



cj : 
cr. 

A wide-spread expectation from a full theory of quantum gravity is the possibility to fix the coefficients appearing 
in the conventional non-renormalizable perturbative expansion seen as an effective field theory (EFT). To address this 
^ • question, a necessary tool is to control the perturbative expansion of the full theory, fn this paper, we investigate 
' this issue in the spinfoam formalism, using the 3d toy model with a single dynamical variable introduced in [ij and 
00 ■ developed in . 

Pursuing a matching with the EFT, while right at the root of many approaches to quantum gravity, most notably 
string theory and the asymptotic safety scenario, has long been obstructed in the spinfoam formalism. This is due to 
CNj ' the difficulty in consistently inserting a background metric to perform the perturbative expansion. The key idea is 
I to relate the n-point functions to the field propagation kernel, via the introduction of a suitable boundary state 0]. 
The boundary state can then be taken to be a coherent state^ peaked on a classical geometry We then expect 
the boundary geometry to effectively induce a semi-classical background structure in the bulk, which allows to define 
the graviton propagator from background-independent correlation functions. 

The structure of this framework is particularly clear in 3d. Considering for simplicity the Riemannian case, the 
spinfoam amplitude for a single tetrahedron is the 6j-symbol of the Ponzano-Regge model, fts large spin asymptotics 



d • is dominated by exponentials of the Regge action for 3d general relativity. This is a key result, since the quantization 
of the Regge action is known to reproduce the correct free graviton propagator around flat spacetime [6|. The role 
of the boundary state is to induce the flat background and to gauge-fix the propagator Q. Thus the framework 
provides a clear bridge to Regge calculus as an effective description of spinfoam gravity. However, there is more to it. 
Indeed, if one works with quantum Regge calculus alone, there are technical problems to go beyond the free theory 
approximation. These are related to the lack of a unique measure for the path integral compatible with the triangle 
inequalities conditions ensuring that the metric is positive definite. The issue is solved in the spinfoam formalism, 
where the triangle conditions are automatically imposed on the 6j-symbol by the recoupling theory of SU(2) and the 
measure is selected by the topological symmetry of the system. Thus the spinfoam approach does reduce to quantum 
Regge calculus at leading order but improves it beyond ^. 
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^ Similar ideas on the use of coherent states lie also behind the study the semiclassical limit in the canonical loop gravity framework Q . 

^ The situation is more complicated in 4d. Developments of this idea have led to the remarkable result that the Barrett-Crane model 
in 4d Riemannian spacetime does reproduce at large scales the scaling behavior of the free graviton propagator (or 2-point function) 
fsl I9I. [Tol llll. [T^. This is crucial evidence towards the correctness of the semiclassical limit of LQG. However the same developments 
also pointed out [sl. Il2ll that the Barrett-Crane model does not reproduce the right tensorial structure of the propagator, thus the model 
fails to reproduce General Relativity in the large scale limit. These results have confirmed the validity of the method, and spurred new 
efforts towards a better understanding of the spinfoam dynamics [T3I. Il4ll . This better behaved models should have a semiclassical limit 
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In this paper we consider the simplest possible setting given by the 3d toy model introduced in [l|, y] and study 
analytically the full perturbative expansion of the 3d graviton. Our results are based on a reformulation of the 6j- 
symbol and the graviton propagator as group integrals and the saddle point analysis of these integrals. We compute 
explicitly the leading order then both next-to-leading and next-to-next analytically and we support these results with 
numerical data. Moreover, it was shown in [3| that deviations of the 6j-symbol from the leading order Ponzano-Regge 
asymptotics do not contribute to the next-to-leading order of the graviton, but that they enter the next-to-next order 
corrections. Here, the exact representation of the graviton propagator as a group integral naturally incorporates these 
deviations. Finally, an interesting side- product of our calculations is a formula for the next-to- leading order of the 
famous Ponzano-Regge asymptotics of the 6j-symbol in the special isosceles configuration. 

In spite of the simplicity of the model, the framework we develop here has rather generic features useful for 
computing graviton correlation functions in non-perturbative quantum gravity from spinfoam amplitudes, although 
it does not allow us to tackle the more general issue of the existence of a relevant boundary state and of the resulting 
EFT-like expansion of the correlation functions for a generic spinfoam triangulation. We nevertheless show that the 
full perturbative expansion of the two-point function in the spinfoam quantization of 3d gravity is computable. We 
hope to apply these same methods and tools to 4d spinfoam models and allow a more thorough study of the full 
non-perturbative spinfoam graviton propagator and correlations in 4d quantum gravity. 

I. THE KERNEL AND THE PROPAGATOR AS GROUP INTEGRALS 

A. The boundary states and the kernel 

Let us consider a triangulation consisting of a single tetrahedron. To define transition amplitudes in a background 
independent context for a certain region of spacetime, the main idea is to perform a perturbative expansion with respect 
to the geometry of the boundary. This classical geometry acts as a background for the perturbative expansion. To do 
so we have to specify the values of the intrinsic and extrinsic curvatures of such a boundary, that is the edge lengths 
and the dihedral angles for a single tetrahedron in spinfoam variables. Following the framework set in [l|, we restrict 
attention to a situation in which the lengths of four edges have been measured, so that their values are fixed, say to a 
unique value jt -\- \ - These constitute the time-like boundary and we are then interested in the correlations of length 
fluctuations between the two remaining and opposite edges which are the initial and final spatial slices (see figure 
[1]). This setting is referred to as the time-gauge setting. The two opposite edges ei and 62 have respectively lengths 
ji + \ and j2 + 5- In the spinfoam formalism, and in agreement with 3d LQG, lengths are quantized so that jt^ ji 
and j2 are half- integers. 




FIG. 1: Physical setting to compute the 2-point function. The two edges whose correlations of length fluctuations will be 
computed are in fat lines, and have length ji + i and j2 + ^ ■ These data are encoded in the boundary state of the tetrahedron. 
In the time-gauge setting, the four bulk edges have imposed lengths jt + | interpreted as the proper time of a particle 
propagating along one of these edges. Equivalently, the time between two planes containing ei and 62 has been measured to 
ber=0i + i)/\/2. 

The lengths and the dihedral angles are conjugated variables with regards to the boundary geometry, and have 
to satisfy the classical equations of motion. Here, it simply means that they must have admissible values to form a 
genuine flat tetrahedron. Note that the dimension of the SU(2)-representation of spin j, dj = 2j + 1 is twice the edge 



given by a modified Regge calculus where the fundamental variables are area and angles, as the one investigated in |15| . 
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length. Setting ke — 4t~^ e — Si,e2, the dihedral angles di, 'i?2 and -dt can be expressed in terms of the lengths 



2d„ 

^,. = 2arccos( /^-^ ) and = arccos f ^=^^=) 



(1) 

provided < 1, a condition ensured by the triangle inequalities. Notice the relation : cos i?t = — cos(^) cos(^). 

We then need to assign a quantum state to the boundary, peaked on the classical geometry of the tetrahedron. 
Since jt is fixed, we only need such a state for ei, peaked on the length ji + 5, and for 62, peaked on j2 + ^- The 
previous works used a Gaussian ansatz for such states. However, it is more convenient to choose states which admit 
a well-defined Fourier transform on SU(2). In this perspective, the dihedral angles of the tetrahedron are interpreted 
as the class angles of SU(2) elements. As proposed in the Gaussian ansatz can be replaced for the edges ei and 
62 by the following Bessel state: 



with 7^ = dj^(l-fc2) (3) 



where A'' is a normalization coefficient depending on 7e. The functions In{z) are modified Bessel functions of the 
first kind, defined by : In{z) — ^ Jq d,(p cos{nip) , and — "i^e/S is half the dihedral angle. The asymptotics 

reproduce the Gaussian behavior peaked around je, with 7e as the squared width: 



1 / 4 u-3.f 

*e(j) = ^y— e -r^ cos{djae). (4) 



The role of the cosine in ^ is to peak the variable dual to j, i.e. the dihedral angle, on the value «£• Then the 
boundary state admits a well-defined Fourier transform, which is a Gaussian on the group SU(2). We parameterize 
SU(2) group elements as 

g{(l),h) = cosipl + i sin (ph ■ a, (/ie[0,27r[, n £ , 

where Ui are the Pauli matrices, satisfying af — 1, and (j> is the class angle of g. Since the group element g{(f),n) is 
identified to g{<p + tt, — n), we can restrict (j) to live in [0, tt). The Fourier transform of ^ is then given by: 

?7=±1 

with = —^—-r sin (dj U + r] ae)) 6"^= ''"'('^+'' (6) 

N sm((p) 

This state is a class function on SU(2), but ^''^''^ alone is not (due to the cj) ^ —0 symmetry refiecting that a SU(2) 
group element and its inverse are simply related by conjugation). The semiclassical analysis is crystal-clear : it is 
peaked around the angle on tt — ag, according to the sign of rj. The sine shifts the mean length to je + \ - 

These states carry the information about the boundary geometry necessary to induce a perturbative expansion 
around it. More precisely, we are interested in the following correlator, 

with OjA3')^^{d^r~dl) (8) 

where the normalisation factor M is given by the same sum, without the observable insertions Oj^ {j'e)- ^^1122 measures 
the correlations between length fluctuations for the edges ei and 62 of the tetrahedron, and it can be interpreted as 
the 2-point function for gravity contracted along the directions of ei and 62. 

The 6j-symbol, as it enters ([7]), emerges from the usual spinfoam models for 3d gravity as the amplitude for a 
single tetrahedron. In the previous work Q we studied the perturbative expansion using its well-known (leading 
order) asymptotics in term of the discrete Regge action (for the tetrahedron). Here instead we use the fact that the 
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6j-symbol for the isosceles configuration admits an exact expression as group integrals,'^ 

■J-^ ? o- r " / dgidg2 X0A9192) Xh{9i92^) X3i{gi) Xh{92) (9) 

J2 Jt Jt) JsU(2)2 

where Xji9) — "in ^^"^ SU(2) character. Then, selecting a specific boundary state as described below, we 
are able to rewrite also ^ as an integral over SU(2). This allows us to study the perturbative expansion as the 
saddle point (or stationary phase) approximation of the integral for large lengths, dj^ » 1. With respect to 0, this 
procedure has the advantage of including the higher order corrections coming from both the Regge action and the 
corrections to the {6j} asymptotics. We will come back to this point below. 

Let us begin by looking at the saddle points of the isosceles 6j-symbol, as the computation of the propagator will 
have a similar structure. We first need the angle of the group elements (7152 and gig^^ : 

"^12 = arccos ( cos ipi cos 02 T " sin 0i sin ^2) (10) 

where we used the notation u — ni ■ 712- Then, expanding the rapidly oscillatory phases in exponential form yields: 

J2 ]\ ^'S} " i ^ I '^'^I'^'^^d^^ /(0i,</'2,^^) e^'-*t»> (11) 

£i,e2,ef2''=l2 = ±l 

with 

f{^^,h,u)^ Bin(0i)sin(02) (^2^ 
sm(0+)sin(0i2) 

*{£}(01>2,U) = (e]*2012 +er2'/'r2) + 2fciei + 2^262 (/>2 (13) 

Let us proceed to the search for the stationary points of the phase ${£}• The variable u only enters (/)f2, and the 
related equation, e]*2^«0i2 + ^i2'^u(f>i2 — Oi is solved hy u — ni ■ n2 — and e'^2 — ^12 ^ ^12- The variational equations 
with respect to 0i and 02 are: 



c^jt(e^2^0.<^i2 + ei2^0e0i2) +c;>ee =0 e = l,2 (14) 



and are solved in [0,7r[ by 



£12 £2 arccos ( ) + (1 + £12 £2)1 

(15) 



-ei2 ei arccos ( z^^^^^ ) + (1 + <^12 ei)| 



fci 

Notice that this result allows us to give a geometrical interpretation to the class angles entering (fTTj) . which is similar 
to the one for the usual integral formula of the squared 6j-symbol: the isosceles 6j-symbol is peaked on half the 
internal, or external, dihedral angles of the classical geometry. Indeed, for example, when £1262 = ~1: the stationary 
angle 0i is 0i = ai = '?i/2, while for £1262 — 1, we have 0i = tt — ai. ^ 

We perform the complete expansion of the isosceles {6j}, using this stationary phase analysis, below in section [TVl 
Now we turn to the graviton propagator. 

B. The propagator as group integrals 

If we insert the expression © into ([7]), the sums over j[ and J2 give the SU(2) Fourier transform of the boundary 
states. Let us first look at the normalization J\f : 

AA = J dgidg2 Xjt (9192^) Xj, {9192) [JI ^-^i (jOXj; (.9i)] [^I *'^2(j2)Xj^(.92)] (16) 
= j dgidg2 Xjt{9i92^) Xjt{9i92) $ei(5i) ^62(^2) (17) 



^ For a general configuration, one has to consider the squared 6j-symbol to have an integral expression. 
This geometric interpretation of the saddle point is not surprising since the 6j-symbol is indeed the unique physical quantum state for a 
trivial topology and a triangulation made of a single tetrahedron [17[ . It satisfies the quantum flatness constraint and can serve as the 
boundary state in the general boundary framework. 
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For large spins, we are interested in evaluating J\f with a saddle point approximation. To that end, let us expand the 
previous expression in exponential form : 



r;i,?)2 = ±l 

with S[^^^y{(f)i,(j>2,u) = i{e'^2't'l2 + ^12'i>12) + E '^'ikef^e{(f'e + Veae) - - kj) Sin^{(l)e + Veae) (19) 

e=l,2 

where the label {e, ry} refers to the dependence on the sign variables, and / is given by (|12p . The crucial point is that 
the phase of (the imaginary part of S^^ rj}) is precisely ^{e} hi p^ . up to constant ai,2 terms which play no 
role in the stationary phase approximation. This means in particular that the imaginary part of S^g^fjy has the same 
saddle points of the isosceles 6j-symbol. 

The same analysis can be performed for the numerator of Wii22- To take into account the observables Oj^, notice 
that the SU(2) character is an eigenfunction for the Laplacian on the sphere with the Casimir as eigenvalue: 

^s^XjW = -4- Xj) = ~{d' ~ 1) X,(0)- (20) 

sm 

This allows to perform the sums over j[ and J2 in ([7|), introducing the Fourier transforms . Again expanding the 
result of these operations into exponential form, one ends up with: 

W1122 = , ^'t'q E ^1 ^2 e^^ f dcjy^dc^^du /(</)i,02,u) e''-^f->(^-^-") 

SZtt^ AI 4 cos-^ Vt _ J 

rn,-n2=±i (21) 

X (^a{,,^}(0i, (1)2) + ^ 
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The functions Oj,;^^}, fe{£,r;} and c^^ ^i} stand for the observable insertions, and are given by 



a 



,n}{4>l,h)^ W { , " SIV? 2((/)e + ??e ae) - 2ke siu 2(0e + ??e Qe)) (22) 



e=l,2 2fc. 



&{e,r,}(0i, 02) = -— cos2(02 + r/2 02) (^-^^ siu^ 2(0i + 771 cki) - 2iei sin2((^i + 771 ai)j + (ei <-> 62) (23) 
C{e,»?}(0i, 02) = 7^ cos 2(01 + Tyi ai) cos 2(02 + ?/2 02) (24) 

fclfc2 

We are now ready to study the large spin expansion of (|2ip . A common choice in the literature is to do so using 
a power series in \/jt (keeping ji/jt and j2/jt constant). However it is more convenient to take as parameter of the 
expansion the dimension dj^ (again keeping ki, k2 fixed). This is the natural choice, as we compute correlations with 
respect to the background geometry with lengths defined by the half-dimensions djj2, and dj^/2. Furthermore, 

as we show below (see (|49|) and ([ST]) ). these are the values of the lengths emerging in the asymptotics of the 6j-symbol. 

As written in pT|) . W1122 corresponds to the mean value of the function a^e H — + for the non- linear 

theory defined by the action Ss^^^rj} ^^nd the integration measure /. The strategy is thus clear: we will compute 
separately the normalisation M and the numerator, perturbatively, with an expansion around the saddle points of the 
action S^^^r)}- 

As stated above, the imaginary part of >5'{e,,,} has the same saddle points of the isosceles 6j-symbol, namely u — Q 
and 0e given in (jisp (independently of 771 and 772)- The extremization with respect to the real part of ^{e^^}, on the 
other hand, constrains the 771 and 7^2 signs. Indeed, for a given solution (0i,02) from ()15p . characterized by £12, ci 
and £2, the signs 771 and 772 have to satisfy: 

sin2(0e + ?7e ae) = 0, fore = 1,2 (25) 

These equations are solved by taking 771 — €2^12 and 772 = £i£i2- This leads to four possibilities, summarized in the 
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following table, 







m = l 


m = -1 


4>i = ai, and 02 = a2, 
^1 — £2 ~ — £12 


01 = TT — ai, and 02 = ^2, 
ei = —£2 — — £12 




(f)i = «!, and 02 = TT — q;2, 
— ei = £2 — £12 


01 = TT — Qfl, and 02 — TT ~ a2, 

£l = £2 = £l2 



The condition e 



12 



"12 



ei2 and rji — £2£i2 and 7]2 = £i£i2 allows us to perforin three sums in (|2ip. the 



configurations for which there is no saddle point being exponentially suppressed : 
1 kik2 



327r2 TV 4 cos2 



E 



£l £2 



ei ,«2,ei2 = ±l 



6id02du /(0i,02,m) e''«'^{=>('^i''^='") 



/ / , , X , ^{e}(01,02) , C{,}(0i,02)n 

(a{.}(0i,02) + + j 



(26) 



and the same for Af without the insertion of ^^^^i^^ ('^{e} + ^{nj/^jt + '^{e}/'^it)' Here the label {£} simply indicates 
the dependence of the functions on the signs £i, £2 and £12- 



II. THE COMPLETE PERTURBATIVE EXPANSION 



The perturbative expansion of the two-point function W is formulated as an asymptotic power series expansion in 
of the type: 



W 



1 



-Wi 



1 



W2 



(27) 



Let us remind the reader that the dimension dj^ defines the length scale of the tetrahedron L = with the 

Planck length Lp = KG. Such an expansion thus matches the typical expansion of quantum field theory correlations 
with quantum corrections ordered in increasing powers of h (and of the coupling constant G) and with L corresponding 
to the renormalization scale. We can thus call the coefficients wi,W2t- the one-loop and two-loop (and so on) 
corrections. 

This perturbative expansion is obtained studying the power series expansion in around each of the four saddle 
points of both the denominator J /expdjjS'^g} and the numerator /(ojj} + Cjjj/d^J/expdjjS'^g}. More 

precisely, we expand the action Ss^^y around its saddle point: the evaluation of 5' at the stationary point gives a 
numerical factor, there is no linear term obviously, the quadratic term defines the Hessian matrix ^{e} and finally 
all the remaining higher order terms (cubic onwards) are kept together to define the potential ^{e}- This potential 
thus contains all higher order corrections to the quadratic approximation to the action S. As such, it does not enter 
the leading order of the two-point function but largely enters its NLO, NNLO and so on, (the loop corrections) as in 
quantum field theory. Then each term in the power series is evaluated as the Gaussian moment with respect to the 
Hessian matrix ^{e} of terms coming from the expansion of / exp dj^rij^-j in powers of dj^. In general many terms 
actually contribute to the same overall order in 1 /dj^ . This intricate structure comes about precisely as in Q because 
the expansion of exjpdj^VL^^i^ gives increasing powers of dj^ while the Gaussian moments have increasing powers in 

On the other hand, a simplification of our calculations comes from the fact that each saddle point gives the same 
contribution. This is a consequence of the symmetry properties of the functions involved under the transformation 
of 0e into vr — 0e. Further, for a given saddle point, the two possible configurations of signs are simply related by 
complex conjugation. The actual sum then ensures the reality of the result. Without loss of generality, we can thus 
restrict the computation to the saddle point (ai,Q;2,0) with ei2 = 1 = — £1 = —£2- There we have : 



a(0i,02)= W (^^^ sin^ 2(0e - ae) -t- 2isin2(0e - Qfe) 



&(01,02) 
C(01,02) 



e=l,2 
1 

1 

fcifc2 



COS 2(02 - a2) 



l-kl . 



2fci 



sin 2(01 — ai) + 2i sin2(0i — ai 



)) + (ei ^ 62) 



cos 2(01 — ai) COS 2(02 — 02) 



(28) 

(29) 
(30) 
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and the potential is extracted from the derivatives of S greater than three, with S given by with the chosen 
signs, 

S((f>i, (l>2, u) — *('?^'i2 + 'Pu) ^ (l ^ ^e) sin^((/'e " cte) + hncar terms (31) 

e=l,2 

Expanding around the background, the inverse of the Hessian matrix is (see the Appendix A for details): 
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l — kf felK2 

feifea 

V ^W^)) 



(32) 



Introducing the shorthand notation 



^p" - (33) 

all possible pairings 

of (ft,...,/32N) 

for /3 € {1, 2, 3}^^, the complete perturbative expansion of the propagator can be written as 



kik2 2d,, 2.^,j=i,2 0'»,a +2.p>2 dP 

In the numerator, the first term gives the leading order contribution in l/dj^ (see next section). It comes entirely 
from the a term in (|28[) . In fact, a and b vanish at the saddle point, and so does the gradient of a, so the expansion 
of a, b and c is dominated by the quadratic term of a. 

All the higher order corrections have been collected in the summations. The coefhcients Np and Wp correspond 
to finite sums: 



^ (2(F + n))!n! V /3 ; f3 J \^,=a,,^^=a^,u=o ^ ' 



and 



"=0 /3e{l,2,3}2<P+") 

^n!l ^ 2 F + 71 ! V Pa ^ ■' ^ Pa J 

n>0 /3';e{l,2,3}2(P+") 

^ (2(P + n- 1 ! V /3b ^•''0bJ 

/3i,e{l,2,3}2(-P+"-i) 



^ 2 P + n-2))! V /3c ^ /3c y J |0i=ai,02=a2,n=O 

/3'cG{l,2,3}2<-P+"-2) 

(36) 



The three lines of p6p are the three separate contributions of the insertions a, b and c. The sum over n defining Wp 
is finite for each of these contributions : n is bounded by 2P — 2, 2P — 3 and 2P — 4 for a, & and c respectively. The 
derivatives of highest order of fl involved in Wp are respectively the 2P-th derivatives, the (2P — l)-th ones and the 
(2P — 2)-th ones, and for Np, the 2(P + l)-th derivatives, all corresponding to n — 1. 

The intricacy of the formulas was anticipated at the beginning of the section. However the reader should be 
reassured that they are simple, if tedious, algebraic expressions. 

The real part 3? in (|35p is consistent with the reality of the initial expression ((2T|) . and arises from the summation 
over the e sign. 



III. THE LEADING ORDER, ONE-LOOP AND TWO-LOOP CORRECTIONS 

We now use psp and ([M)) to obtain explicitly the first orders of the expansion. The leading order (LO) and the 
next to leading order (NLO) have already been obtained in a quite different way in We here recover them quickly. 
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The computation of the next to next to leadmg order (NNLO) is then completely new. It is shown in Q that the 
NNLO needs the corrections to the asymptotics of the 6j-symbol, i.e. to the Ponzano-Regge formula. The success of 
our method resides in the fact that such corrections are naturally contained in the exact group integral expression ^ 
of the kernel. 

The LO is obtained evaluating the normalization at the saddle point, /o — /(ai,Q;2,0), and the numerator at the 
non-zero second derivatives of a: 



"1122 — 



sin(2dj, - l)§t 



dj^ cos-dt sin(2djj - 



(37) 
(38) 



This reproduces the expected l/dj^ scaling behavior of the LO. The difference in the coefficient with [l|, Q comes 
from the different boundary state used. In particular notice that while 'dt{_ki, is a constant, the dependence upon 
of the second fraction produces spurious oscillations. These can be reabsorbed in the boundary state, replacing 
sva.dj^{4) + 77 tte) in ^'e with sin {dj^{(f) + rj ae) + dj^dt)- The Fourier transform is then: 



«'e(j) 



N 



b-ie|(y) cos(djae + dj^-dt) 



d,Jt) 



(39) 



This does not affect the asymptotic behavior of ^'e(i)- With this replacement, we obtain the same result of (cf. 
equation (37)) for the isosceles case. 



"1122 



1 sm^-dt 
dj^ cos dt sin^t^t 



(40) 



Even if the LO now matches the previous results presented in [3| , the higher orders will differ because of the different 
boimdary state used. 

For the sake of a simpler presentation, we will report the NLO and the NNLO for the equilateral tetrahedron, 
ki — k2 — 1/2 and — arccos— ^. The general expressions in terms of fci and k2 are indeed quite cumbersome. This 
choice will also facilitate the comparison with numerical simulations of ([7]). 

The NLO is then obtained from the coefficients Ni and W2. To keep compact expressions, we adopt the following 
symbolic notation for the contractions of derivatives with Gaussian moments : for functions / and h of 4>2 and 
u, define: /„ h„, A'l^ = Ezi,... ,z„=i,2,3 Eji,- j„=i,2,3 cfi,.. ^n,-,,J^ ^(^i,- ji,- j™) evaluated at the 
saddle point (ai,a2,0) with ei2 — — ci = —£2 = 1- For example, the LO of (|34p can be written ^02^^^. In A^i, 
three powers of f2 appear, 17°, 17 and il^. Using the boundary state (|39p. we have: 



We proceed in the same way for the three contributions to W2: 



h A^^ + (/i ^3 + /o Si)A-^^ + ^^3^3 A-^^ 



fo 



(02 /2 + 03 /l + /o 04)^44 ^ + (02 fl S3 + /o(a3 5*3 + a2 Si))AQ ^ + y 02 5*3 53 ^ 

+ {fo 62 + 61 /i)^2"' + fo bi S3 A^' + fo Co 



(41) 



(42) 



After straighforward algebra we obtain the NLO, of order l/d^ 



wl 



NLO 
122 



1 



511 

432 4 



(43) 



These results for the LO and NLO are well-confirmed by numerical simulations, as one can see from figure [21 An 
agreement with 0.58% of error for the LO, and with 1.7% error for the NLO is reached between the coefficients of 



these orders for dj^ — 201 (i.e the representation of spin jt — 100). 



All orders of the expansion can be computed using the above recipe. From this point of view, the NNLO is of no 
particular specificity. We need the expansion of the action (or equivalently fl) until the sixth order. The highest order 



FIG. 2: Log-log plots comparing numerical simulations with analytical results. Left plot: a numerical simulation of ^ 
(diamond symbol) compared with the leading order of (|43|) (star symbol). Middle plot: the next to leading order of (|43p . in 
star shape, with the numerics, in diamond shape. Right plot: the next to next to leading order (|48|l . 



correlator A^^ is of order 12 for the normalisation, and respectively 14, 10 and 6 for the insertion of a, b and c. 



5r( 



h ^4~' + (/3 ^3 + f2Si + h S5 + fo Se)A^^ + -{h S3 S3 + 2/1 ^3 ^4 + fo ^4 ^4 + 2/o ^3 S5)A^^ 



3! 



(/l ^3 ^3 S3 + 3fo S3 S3 54) A 
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fo 
4! 



S3 S3 S3 S3 A-^2 



(44) 



We also write Wg = ^^3^"^ + W^^''^ + W^''\ with: 



W. 



{fo ae + /i 05 + /2 0,4 + fi as + /4 02) A^^ + {{fo + /i a4 + /2 as + /s a2)S3 

+ {fo 04 + fl 03 + /2 a2)'5'4 + (/o 03 + /i 02)55 + fo 02 56)^^^ + ^((/o 04 + /l 03 + /2 02)^3 ^3 

+2(/o 03 + fl 02)53 54 + /o 02(54 54 + 253 S5))A^^ + ^((/o 03 + /i 02)53 53 53 

+ 3/0 0,2 S3 S3 54)Aj^2^ + ^"^2 53 53 53 53 A^l 



w. 



(b) 



{fo hi + /i &3 + /2 62 + /3 hi)Al^ + ((/o 63 + /i 62 + /2 61)53 + (/o h2 + /i 61)54 
+/o 61 55)^6-1 + i((/o 62 + /i 61)53 53 + 2/0 61 53 54)^3-^ + ^61 53 53 53 A~^ 
14^3^^) = 5R(z e^'^' [(/o C2 + /i ci + ^ co)^^-! + {{fo ci + /i co)53 + fo co 54)^-1 + ^cq 53 53 ^J) 



The NNLO is thus computed to be: 



NNLO 
1122 



511 



520507 



d,-, 432 S 157464 £ 



(45) 
(46) 

(47) 
(48) 



This result is again supported by numerical simulations, see figure [2] An agreement with 11.3% of error is obtained 
for the coefficient of l/dj^ at dj^ — 201. The error can be reduced pushing the simulations to higher values of dj^. 

IV. PERTURBATIVE EXPANSION OF THE ISOSCELES 6J-SYMBOL 



The procedure described above can be applied directly to the isosceles 6j-symbol (j9l), obtainingthe known Ponzano- 
Regge formula and its corrections. This is interesting for a number of reasons. As discussed in [4I, the corrections to 
the Ponzano-Regge formula are a key difference between the spinfoam perturbative expansion studied here, and the 
one that would arise from quantum Regge calculus. The 6j-symbol is also the physical boundary state of 3d gravity 
for a trivial topology and a one-tetrahedron triangulation. In 4d, it appears as a building block for the spin-foams 
amplitudes, such as the 15j-symbol. Thus, with regards to many aspects of spin-foams in 3d and 4d, in particular for 
the quantum corrections to the semiclassical limits, it would be good to have a better understanding of this object 
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beyond the Ponzano-Regge asymptotics. This is what we do in this section, performing a perturbative expansion of 
the exact expression ^ for the isosceles 6j-symbol. Indeed, it is a simpler application of the procedure developed 
above for the propagator. 

As written in (|lip . the 6j-symbol is the partition function for the theory defined by the action ^{e} and the 
integration measure /. Notice first of all that the Regge action of 3d gravity Sr — ^ dj^^ naturally appears as the 
evaluation of ^{e} at the saddle points: 

" " (49) 
SR = d,^ {2^t + 2kiai + 2k2a2) 

We then proceed exactly as for the propagator, knowing that for each configuration of signs, the {6j} is peaked on 
the classical geometry of the tetrahedron. The perturbative expansion with respect to this flat geometry is thus given 
by the Gaussian moments of the Hessian matrix -ff{e} of Let us stress that, in contrast with the previous studies 

of the asymptotics of the 6j-symbol, we are here scaling the lengths of the tetrahedron (or equivalently dj^), keeping 
the length ratios ki and k2 fixed, instead of scaling jt. 

As for the propagator, the four saddle points give the same contribution, and the two sign configurations of a 
given saddle point are related by complex conjugation. This can be done in a quite explicit way. Introduce lo to 
be the truncated Taylor expansion of starting at order three onwards, around the saddle point (ai,Q;2,0) with 
£12 — —ei = —62 = 1- Let H^^ be the corresponding inverse of H^^y. 



2 



/ r-^COt^t , , / , \ 

1 Vcoti?t 



(50) 



^l-fc^y'l-fcl siniJt l-fe; 

V " " i-ikf+k-i)/ 

We also introduce the volume of the tetrahedron, which enters the Gaussian integrals of H: 

yt = ^ k^k2^i-{ki + ki). (51) 

The expansion of this isosceles 6j-symbol is then (see appendixlBl for more details): 
where the coefficients Cp, for P = 2p, 2p + 1, are given by finite sums: 

n=0 /3e{l,2,3}2(-P+'«) 

Thus, all even orders are in phase with the leading order asymptotics, given by the original Ponzano-Regge for- 
mula in cos{Sr + 7r/4). This leading order is easily recovered by computing the coefficient Cq, with /(ai,Q;2,0) = 

yr^yr^: 



^ cos (Sr + ^) . (54) 



J Ji Jt I ^ 

\j2 Jt jt J LO ^/viWt ""^ ' 4 

On the other hand, all odd orders are out of phase (or in quadrature of phase) with this leading order. If we were 
scaling the spin jt instead of the length djt/2, the result would not have had such a simple structure with sines and 
cosines being mixed up at all orders (but leading). 

This asymptotic series formula for the isosceles tetrahedron shows that only the Regge action is relevant and no 
other frequency appears in the 6j-symbol. We believe this feature to generalize to the generic 6j-symbol since its 
asymptotics can also be extracted using saddle point techniques [IQ] . 

The coefficient of a given order is simply given by the contractions of the derivatives of /a;" with the Gaussian 
moments. For a given order P, the highest derivatives of uj involved correspond to n = 1 in (|53p and equals 2(P+ 1). 
For instance, the NLO is obtained by setting P = 1. With the notations of the previous section, we have: 

Cl = /2 - ifl LO3 + fo L04)H^' + LU3 H,' (55) 
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FIG. 3: Differences between the 6j-symbol and the analytical result (|56|l for three pairs (^1,^2) — {k,k): from left to right, 
k — 1/2, k = 3/14, k = 3/42. The X axis stands for dj^ = Ndj^, for djg respectively fixed to 1, 7, and 21, while A'' goes from 
200 to 800. 



and introducing the reduced volume v — Vt/d^^ : 



Ji Jt Jt 



1 / ""X cos^i^t Pi{ki,k2) 

cos ' • 



[j2 jt Jt] NLO ^mVt V 4/ dj.^nWt 48(12t-)3 
where P(fci, fc2) is a symmetric polynomial in fc^ and 
Pi(fci,/c2) =3(l-fc?)2(i-2fc2) + 3(l-fc2)2(i_2fc2)--3 + 46fc?fc^-t-25fc;'fc|--44(fc^ (57) 

This polynomial is not simply related to the volume and we haven't succeeded in providing it with a geometric 
interpretation. It would nevertheless be very interesting to understand its geometrical origin in order to interpret 
physically the higher order corrections to the graviton propagator. 

For extremal values of A;i, this polynomial simplifies. We get ^1(0, k) = 3(1 — ^2)^(1 — 2k'^) for fci = 0. At the other 
end at ki = 1, we obtain Pi(l,/c) = — 4A:'^(1 — fc^) with obvious roots and 1. Let us point out that ki^2 actually 
never physically reaches these extreme values and 1, but its bounds depend on the representation jt (due to the 
SU(2) recoupling theory): 

^— < fee < 1 

2d,, - - 2d,, 

When ki reaches these extreme values, the coefficients of Pi are polynomials in l/d|j- ^ 

The result (j56p is confirmed by numerical simulations, see figure [31 These plots represent numerical simulations of 
the {6j} minus the analytical formula ([56)l . for three pairs (/ci, k2)- We have used in these simulations the particular 
case ki = k2 — k, for which 

Pi(fc,/c) = {l-k^)(3-21k^ + 55k-^ -A5k^) (58) 



whose only root in [0,1] is fc = 15[(10(81\/3T0 + 1450))3 + (10(81^310- 1450))3 -1-55] « 0.8248 thus inducing 

the vanishing of the NLO. To enhance the comparison, we have multiplied by dj, to see how the coefficient of the 
NLO is approached, and suppressed the oscillations by dividing by those of the NNLO, cos{Sr + j). The numerics 
support both the coefficient and the phase. 

Notice that in the equilateral situation, fci = fc2 = 1/2 represent (half) the length ratios but also (half) the spin 



For instance, when the edge ei is at minimal length, ji = or fci = 7^ — , the coefficients of Pi read: 

, 1 , , 1 5 1 X , 23 11 5,0, 11 25 , 4 , 5 , 

Pi( ,fc2 = 3(1 - ^ + ; j^) + ( - 12-1- ^= i- + f^)fc? -I- (15 - ^ 4- 7r)kt+ ( - & + 

2di/ ' ^ S. led-* 32d& ' ^ 2d2 32^6 > ^ S. led* > ^ ^ 2d^ ■ 

Jt Jt Jt Jt Jt Jt Jt Jt Jt 
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ratios. We can thus switch easily to the usual \/jt expansion: 



NLO 



25/4 



-4 



: cos 



7^^ 

A' 



31 . /„ TT, 

-== sm [Sr + -) 

72 -2^^^ ^ 



cos {Sr + f ) 



3 cos + -) + sm (5;^, + -) 



(59) 
(60) 



This point of view shows that it is much more natural to study the asymptotics of the 6j-symbol in term of the inverse 
length \/dj instead of the inverse spin label 1/j. For instance, the leading order coefficient is given in term of the 
volume Vt of the tetrahedron with edge lengths given by the dj's and not the fs. 

Finally, we point out that the asymptotics given above in term of the cosine and sine of the Regge action holds for 
mid-range values of fci, k2 and it breaks down for fci, k2 close to their extremal values and 1. Indeed when fc2 = the 
asymptotics are given in term of Airy functions while when ^2 = 1 they are given by the (non-oscillatory) exponential 
of the Regge action. The interested reader can find details and references in [l3|. 



Conclusions 



We have shown it is possible to compute analytically the two-point function - the graviton propagator - at all orders 
in the Planck length for the 3d toy model (the Ponzano-Regge model for a single isoceles tetrahedron) introduced in 
[l[ . This builds on the previous work [2] where the leading order and first quantum corrections were computed using 
the asymptotics of the 6j-symbol in term of the Regge action. Here, we introduced a representation of the relevant 
6j-symbol and of the full graviton propagator as group integrals over SU(2). Then one obtains the expansion of the 
two-point function as a power series in the inverse spin label (or equivalently in the Planck length) by expanding 
these group integrals around their saddle points. We computed explicitly the first and second order corrections to the 
leading order behavior and matched them successfully against numerical simulations. 

A side-product of these calculations is the corrections to the Ponzano-Regge asymptotic formula for the 6j-symbol 
for an isosceles tetrahedron (when four representations are taken equal). We obtain a series alternating cosines and 
sines of the Regge action for the tetrahedron (shifted by 7r/4). We computed explicitly the next-to- leading order 
correction and checked it numerically. An open issue is the geometrical interpretation of the polynomial coefficient 
Pi{ki,k2) in front of this first order correction. 

To conclude, we have shown how to carry out the calculations of the spinfoam graviton propagator at all orders at 
least in this simple setting. We hope to apply the present methods and tools to more refined 3d triangulations [lB| 
and to compute spinfoam correlations for 4d quantum gravity along the lines of [111, [l^, [3, [3 ■ 
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APPENDIX A: DETAILS FOR THE PROPAGATOR EXPANSION 



The key object containing the quadratic fluctuations and their corrections is the generating function Z(J), which 
is the Gaussian integral of the Hessian matrix A with a source J, evaluated at each saddle point: 

(Al) 



with X = {(f)i, (f)2, u). The Hessian matrix is given, for all configurations of signs, by: 



A 



' sin -dt 
2 ^1^: 



1_ ^2ei€2€l2(f -^?t) 2 ^^^^ e~ 

sin COS 'dt 



V 



sin COS 'dt 





sin i?£ 



,ieie2ei2(f -i?t) 









-2^e,e2e,2 [l - (kf + kj)] 



(A2) 
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Due to the initial symmetry between the groups elements gi and g2, clearly expressed in (|17p . ^{e} is invariant under 
the reversing of ei and €2- ^{e} has also the property of transforming into its complex conjugated matrix when 
reverting all signs, an operation which does not change the saddle point (0i, 02, 0) considered, and equivalently under 
the flipping of ei2. A straightforward calculation yields: 



3/2 



with Z. 



f- 



-16162612 



(l-fc2)(l-fc2)^2|cost9t| 



and A^^]^- 



ei62 



V 



, COS ^iei€2€i2^t 





cos i?t „ieie2ei2i5t 
fclfc2 

1 








^tan?9t e*'^!'^'!^^ 



(A3) 
(A4) 

(A5) 



and benefit from the previously mentioned symmetries of ^{e}- The symmetry flipping 61 and 62 means 
that the saddle points (ai,Q;2,0) and (tt — ai,7r — 0^2,0) have the same Hessian matrices, for a fixed 612. Moreover, 
flipping 612 while going to the saddle points (ai, tt — 0:2, 0) and (tt — ai, 0:2, 0) does not change A^^y and Z^^y, up to 
a change of sign for the non-diagonal coefficients of ^{e}- 

Let us focus on the normalization Af. The Gaussian moments are generated by successive derivations of Z(J) with 
respect to the source, and they are contracted with the derivatives of f exp dj^Q^^y, which we expand into powers of 



djj. We thus have: 



327r2 ^ 



eie2 -^{e} e 



djtS{,}(01,02,O) 



E 

/3e{l,2,3}^ 



n> 



(2^)! C 



W^d}^if^h)^^.._$.,A-:yj (A6) 



where the correlators - are defined according to Wick's theorem: 



all possible pairings 
of (/3i,...,/32jv) 



■ "^{«},ft2JV- 



(A7) 



As is a Taylor expansion into powers of {(jji — 0i), {(f)2 — (t>2) and u, whose minimal order is 3, the power n of f2 
in (|A6|) is bounded from above by N: 3n < 2N, and the sum over n is thus finite for each N. The power of 
receives two contributions: one, positive, from the Gaussian moments, and the other, negative from the expansion of 
expdj^fl^^y. We can identify the coefficients of a given order by the simple change of variables P = N — n. Introducing 
the explicit expressions of Z^^y and S'{e}(0i, 02, 0): 



2P 



with 



Up 



E 



-1 



E up 
dP 



32(1 - kl){\ - fci)V2^|cosi9t| d][' 



(A8) 



E E 2^ 46162612 e 

n=0 ^gjl 2,3}2(-P+i.) ei,e2,«l2 



■ieie2ei2(2d3t - 



= 51' 

(2(F + n))! n! /3 



2(P+n) 



{t} ) 101,02 
M=0 



Let us further simplify the coefficients wp by performing the sums over 61, 62 and 612. First, notice that the sign of 
the imaginary part of ^{^y is 6162612. Since / is real, and considering the symmetry properties of | given in (jASp . 
it is clear that when the signs 61, 62 and 612 are all flipped, the derivatives are evaluated at the same saddle point and 
Mp is transformed into its complex conjugate. Thus, let us work with a fixed value of 6162612, say 1, and consider the 
basic properties of the functions / and S^^^t^ minus the linear parts in 0i and 02 (its derivatives greater than three 
are those of f2). More precisely, we are interested in how these functions and their derivatives, evaluated at a given 
saddle point, transform when the saddle point is changed. Let us see for instance the differences when going between 
the saddle points (0i = ai, 02 = Oi2) and (0i = tt — ai, 02 = 012). 

Having impose the value of 6162612, this change of saddle point is determined by the flips of 61 and 612. We have: 
/(7r-0i, 02, u) = /(01, 02, u) and (0^2 + 0r2)(''' - <?^i, <^2, u) = 27r- (0+2 +0r2)(0i, 02, u), while the real part of S{^y is 
non-zero only when derivated an even number of times. Thus, /r2"(0i, 02, u) equals f^l^in — 0i, 02, u) when we flip 
in the same time 612 in front of (0^2 + 0r2) in S^^y. This means that each derivation with respect to 0i flips the sign 
between the two saddle points considered. There is now three possibilities: (i)such a derivation is contracted with 
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another derivation w.r.t. 0i through A^^^ = 



4(l-fc?) 



then the sign is changed twice, i.e. there is no change of sign. 



(ii)It is contracted with a derivation w.r.t. u via A^^ which is zero, so that there is in fact no contribution. (iii)It is 
contracted with a derivation w.r.t. 02 via A^2 whose sign changes under the flip of ei. Thus these two saddle points 
give the same contribution. The proof can be repeated between the four saddle points. Finally: 



J\f- 



-1 



Np 



4(1 - kf){l - fci)V27r|cosi9t| df^' € 



(AlO) 



with Np given by ((35)) . 

The same analysis can be performed for the numerator of the propagator. One has simply to take into account 
the fact that the insertion ^^^^2^^ ('^{e} + ^{e}/'^it + '-{ej/'^jt) involves three different powers of dj^. To perform the 
sums over the signs ei, €2 and ei2, first notice, like for the denominator, that flipping of all them three transforms the 
coefficients into its complex conjugate. Then, restricting attention to a fixed value of the product £162612, it is easy to 
check that the derivative of a, b and c w.r.t. evaluated at (/)i = tt — ai is equal to (—1)^^ times that evaluated at 
= Qfi, while flipping ei and £12, with pi being the number of derivatives w.r.t. (f)i. The same is true for 02. Thus, 
we can reproduce the previous argument showing that the four saddle points give the same contribution. This leads 
us to: 



W^1122 = 



-kik2 



16(l-fc2)(l 



A:|) cos^ y'27r|cos 'dt 



M d 



i3/2 



P>2 "jt 



with Wp given by 



APPENDIX B: DETAILS FOR THE EXPANSION OF THE 6J-SYMBOL 



(All) 



Let us compute the generating function: 



Zj^f(j) 



dX e-^^"i'i^+J^ 



with H 



{^} 



2i£l2 

siwdf 



^/T^^/T^ £162(1 -fci) cos i9t 



\ -6162 [1 - (fc? + fc|)] cos 1?t 

Taking care of the fact that i?{e} has purely imaginary coefflcients, one has: 



with Z 



(6j) 



and H, \ 



«6l2 



1-k 



^l-kl^l-kl sin i?t 



V 







cot dt 




where the volume Vt is given by ([5T|) . Using (H^ . we obtain an expression similar to (|A8[ 



tan?? 



(Bl) 
(B2) 



(B3) 
(B4) 

(B5) 



Iji Jt Jt 



1 



E 



[j2 Jt jt j S^l-k'f^l-k'i^UWt ^0 dP 



with the series coefficients in term of the Hessian: 

2P 



cp=y: 



1 



(2(P + n))! n! 

n = ^ ^ " £1,62,612 



J2 (*^12)" e-''i'^'i2(^« 



+ t) 



E 



/3G{1,2,3}2(-P+") 



^2(P+n) / ^ n \ _ TT-1 



(B6) 



(B7) 
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We are now in position to repeat the arguments of the previous section. The symmetries of the functions /, iei2'-^, 
combined with those of H'^J^ imply that the four saddle points contribute the same. Moreover, the two configurations 
of signs corresponding to a given saddle point, which are related by flipping ei, 62 and ei2, are related by complex 
conjugation. The coefhcient Cp is thus completely determined by the saddle point (ai, a2, 0) with £12 = ~ei — —ti = 
1. Writing B^^^ = «ei2i?{l}, we have that H^^^^ - = (iei2)^+"i?~/j ^, and: 

2P / ,x„4.|£| 

V / ^ (2(P + n))l n' ^ ''|(ai,a2,o) /? I'^i^- ^i- '=2-1 v J 

n=0 /3G{l,2,3}2(f'+") 

It is then clear that 3? (i^^^LfJ e^'^'^n+f ) J is simply cos {Sr + f ) for even P, and sin {Sr + f ) for odd P. 
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